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This talk summarizes results for all the scale independent renormalization constants for bilinear currents {Za, 
Zy, and Zs/Zp), the improvement constants (ca, cv, and ct), the quark mass dependence of Zo, and the 
coefficients of the equation of motion operators for 0{a) improved lattice QCD. Using data at /3 = 6.0, 6.2 and 
6.4 we study the scaling behavior of these quantities and quantify residual discretization errors. 



The use of axial and vector Ward identities has 
proven to be a very efficient and reliable way of 
extracting the improvement and renormalization 
constants for the 0{a) improved fermion action. 
The methodology, references to previous calcula- 
tions, and the notation we use are given in [Q. 
The new features of our calculation summarized 
here are: new determinations of ca including 
0{vn?'a?') corrections and using non-zero momen- 
tum correlators; improved chiral extrapolations 
in the extraction of ct and bp — bA', and a 
quantitative comparison of the scaling behavior 
of the differences between our results and those 
of the ALPHA collaboration ||J^ and 1-loop 
perturbation theory. Results are summarized in 
Table |l| and will be presented in detail in Q . 

The first feature we discuss is the need for 
including an O(m^a^) term in the extrapolation 
of CA to the chiral limit. Data at /3 = 6.4 is 
shown in Fig. Table |l| gives results from both 
the linear and the preferred quadratic fit. Our 
results show a weak dependence of ca on P in 
the range 6.0 — 6.4, unlike that found by the 
ALPHA collaboration, but consistent with the 
recent results by Collins et al. 

The second new feature is the demonstration 
that consistent estimates for ca are obtained 
from correlators with zero and non-zero momen- 
tum once additional 0{p^a^) errors are accounted 
for. A plot of CA versus (12pa/7r)^ at /3 = 6.4 is 
shown in Fig. |[ We find that a linear extrapola- 
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Figure 1. Comparison of linear and quadratic ex- 
trapolation of CA to the chiral limit 

tion to p = Q yields results consistent with those 
obtained using zero momentum correlators, and 
with a slope of expected magnitude. 

The second point concerns the chiral extrapola- 
tion for Z^, Ct, bp — bA- Our estimates presented 
in were based on constant fits as these quanti- 
ties are not expected to have 0{ma) corrections 
if the theory is fully improved to 0{a). We now 
advocate using results of linear extrapolation 
(marked with an asterisk in Table as our data 
show a dependence on m. Such behavior can 
be explained by O^aAgcD rna) corrections which 
can arise as a result of using a mass-dependent 
CA in intermediate stages of the calculations. To 
show the size of this effect, we give both estimates 
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Table 1 



The first error in LANL estimates is statistical, and the second, where present, corresponds to the 
difference between using 2-point and 3-point discretization of the derivative in extraction of ca- Asterisks 
mark values which include 0{ma) corrections in the chiral extrapolations. 





13 = 6.0 


P = 6.2 


/? = 6.4 




LANL 


ALPHA 


P. Th. 


LANL 


ALPHA 


P. Th. 


LANL 


ALPHA 


P. Th. 


csw 


1.769 


1.769 


1.521 


1.614 


1.614 


1.481 


1.526 


1.526 


1.449 




+0.770(1) 


+0.7809(6) 


+0.810 


+0.7874(4) 


+0.7922(4)(9) 


+0.821 


+0.802(1) 


+0.8032(6)(12) 


+0.830 


z\ 


+0.807(2)(8) 


+0.7906(94) 


+0.829 


+0.818(2)(5) 


+0.807(8)(2) 


+0.839 


+0.827(1)(4) 


+0.827(8)(1) 


+0.847 




+0.802(2)(8) 






+0.815(2)(5) 






+0.822(1)(4) 






^0 / -70 


+0.842(5)(1) 


N.A. 


+0.956 


+0.884(3)(1) 


N.A. 


+0.959 


+0.901(2)(5) 


N.A. 


+0.962 


CA 


-0.037(4)(8) 


-0.083(5) 


-0.013 


-0.032(3)(6) 


-0.038(4) 


-0.012 


-0.029(2)(4) 


-0.025(2) 


-0.011 


CA* 


— 0.038(4) 






— 0.033(3) 






— 0.032(3) 








-0.107(17)(4) 


-0.32(7) 


-0.028 


-0.09(2)(1) 


-0.21(7) 


-0.026 


-0.08(1)(2) 


-0.13(5) 


-0.024 


Ct 


+0.063(7)(29) 


N.A. 


+0.020 


+0.051(7)(17) 


N.A. 


+0.019 


+0.041(3)(23) 


N.A. 


+0.018 


CT* 


+0.076(10) 






+0.059(8) 






+0.051(4) 






bv 


+ 1.43(1)(4) 


N.A. 


+ 1.106 


+ 1.30(1)(1) 


N.A. 


+ 1.099 


+ 1.24(1)(1) 


N.A. 


+ 1.093 


bv 


+ 1.52(1) 


+ 1.54(2) 


+ 1.274 


+ 1.42(1) 


+ 1.41(2) 


+ 1.255 


+ 1.39(1) 


+ 1.36(3) 


+ 1.239 


bA - bv 


-0.26(3)(4) 


N.A. 


-0.002 


-0.11(3)(4) 


N.A. 


-0.002 


-0.09(1)(1) 


N.A. 


-0.002 


bA - bv 


-0.24(3)(4) 


N.A. 


-0.002 


-0.11(3)(4) 


N.A. 


-0.002 


-0.08(1)(1) 


N.A. 


-0.002 


bp — bs 


-0.06(4)(3) 


N.A. 


-0.066 


-0.09(2)(1) 


N.A. 


-0.062 


-0.090(10)(1) 


N.A. 


-0.059 


bp - bA 


-0.07(4)(5) 


N.A. 


+0.002 


-0.09(3)(3) 


N.A. 


+0.001 


-0.12(2)(5) 


N.A. 


+0.001 


bp - bA* 


-0.08(30) 






+0.03(10) 






-0.02(4) 






bA 


+ 1.17(4)(8) 


N.A. 


+ 1.104 


+ 1.19(3)(5) 


N.A. 


+ 1.097 


+ 1.16(2)(3) 


N.A. 


+ 1.092 


bA 


+ 1.28(3)(4) 


N.A. 


+ 1.271 


+ 1.32(3)(4) 


N.A. 


+ 1.252 


+ 1.31(2)(1) 


N.A. 


+ 1.237 


bp 


+ 1.10(5)(13) 


N.A. 


+ 1.105 


+ 1.23(11)(7) 


N.A. 


+ 1.099 


+ 1.13(4)(7) 


N.A. 


+ 1.093 


bs 


+ 1.16(6)(11) 


N.A. 


+ 1.172 


+ 1.31(10)(6) 


N.A. 


+ 1.161 


+ 1.22(4)(8) 


N.A. 


+ 1.151 



m Table |1|. Even though a fit linear in ma removes 
only part of the O(a^) corrections, we choose it as 
our preferred value as it is less sensitive to the m 
values used in the fit. One exception is 6p — fe^ at 
f3 — 6.0 for which the constant fit is our preferred 
estimate as the data do not show a linear term. 

Our estimates of Z^, Zy, ca, cv differ signif- 
icantly from those obtained by the ALPHA col- 
laboration who used the Schrodinger functional 
method, while those for by agree. We expect 
the difference (for example AZ^ = Z^{LANL) - 
Z^ {ALP HA)) to vanish as 0{a^) for and Z^, 
and as 0{a) for ca and cy. To test this we fit the 
data assuming this leading behavior and requir- 
ing that the difference vanish at a = 0: 

AZ^ = -(55a)2 -(464a)3, (1) 
AZ^ = -(246a)2 + (629a)^ (2) 



AcA = -(18fa) + (763a)2, (3) 
Acv = -(367a) + (669a)2, (4) 

where a is in units of (MeV~^) and has values 
1/2120, 1/2910 and 1/3850 at the three p. Con- 
sidering that the expected size of the terms is 
0{aAQCD)^ , all the coefficients look reasonable, 
however, the errors in them are large. We make 
the following observations: 

• The difference in Zy is dominated by the 
O(a^) term. 

• The errors in the coefficients for Z^, and cy 
are > 100%. This is not surprising since the 
combined error at each /3 is approximately 
equal to the difference. 

• The coefficients in the fit for Aca have rea- 
sonable errors, however the fit is dominated 
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Table 2 

Results for off-shell mixing coefficients. 



- 




6.0(f) 


6.0(b) 


6.2 


6.4 




C'y + C'p 


2.82(15) 


-^2.68(19) 


2.62(8) 


2.44(4) 




+ 


2.43(24) 


-H2. 12(31) 


2.43(14) 


2.27(6) 




2c'p 


0.88(97) 


-0.65(57) 


1.82(24) 


1.85(8) 




Cg + c'p 


2.44(13) 


+2.40(13) 


2.40(7) 


2.27(4) 






2.40(18) 


-H2. 27(20) 


2.42(9) 


2.28(5) 






2.38(50) 


-1-3.00(37) 


1.72(16) 


1.52(4) 






1.99(56) 


-1-2.45(46) 


1.53(20) 


1.35(6) 






0.44(49) 


-0.33(29) 


0.91(12) 


0.93(4) 




c's 


2.00(48) 


+2.72(33) 


1.49(14) 


1.35(4) 






1.96(49) 


+2.60(38) 


1.51(15) 


1.36(4) 



Figure 2. Evidence for additional O(p^a^) errors 
in ca extracted from p ^ correlators. 

by the significant difference dX (3 — 6.0. We 
note that the (3 dependence of the difference 
is similar for ca, Z^a^ ^^'^ ■ 

Using our three data points we can also fit the 
difference between the non-perturbative and tad- 
pole improved 1-loop estimates as a function of 
the leading residual discretization error in a and 
perturbative, 0{o?^), corrections. The results are 





-(158a)2 


- (1.4a,)2 


(5) 


AZ" = 




- {lAasf 


(6) 


^zliz% = 


-(502a)2 


- (1.8a,)2 


(7) 


ACA = 


-(13a) 


- (1.3a.)2 


(8) 


Acv = 


-(51a) 


-(1.7a.)' 


(9) 


Act = 


(94a) 


+ (0.8a,)2 


(10) 


Abv = 


(930a) 


- (2.6a,)2 


(11) 


Abv = 


(429a) 


+ (1.5a,)2 


(12) 



where a, expressed in MeV~^, — g'^/{A'KUQ) is 
the tadpole improved coupling with values 0.1340, 
0.1255 and 0.1183 at the three /?, and uq is 
{plaquettey^^. The errors in the other b are too 
large to allow any meaningful fits. 

The errors in the coefficients for the three AZ's 
are reasonably small, providing some confidence 
in the fits. Over this range of (3, the perturbative 
corrections dominate the differences in Z^ and 



Zy, whereas in Zp/Zg, the two corrections are 
comparable. 

The errors in the coefficients for the three Ac's 
are large. Even though the coefficients are of the 
size expected, it is important to note that the 
non-perturbative estimates are 2 — 4 times the 
perturbative values. 

Both corrections are large in Aby and Aby, 
with the discretization error being the larger of 
the two. 

Overall, these fits, since they are based on data 
at just three (3 values with 1/a between 2.1 and 
3.86 GeV and since we have ascribed no errors 
to a or as, should be considered indicative and 
qualitative and certainly not sufficient to draw 
precise conclusions. This is why we refrain from 
quoting errors in the fits. 

Finally, in Table |^ we present results for the 
coefficients of the equation of motion operators. 
Estimates aX (3 = 6.0 are poor, but become rea- 
sonably precise at /3 = 6.2 and 6.4. We find that 
except for Cp, the corrections to the tree level 
value c'q — 1 are large. 
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